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The properties of warm dilute alpha-nucleon matter are studied in a variational approach in the 
Thomas-Fermi approximation starting from an effective two-body nucleon-nucleon interaction. The 
equation of state, symmetry energy, incompressibility of the said matter as well as the a fraction 
are in consonance with those evaluated from the virial approach that sets a bench-mark for such 
calculations at low densities. 
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I. INTRODUCTION 

Cold nuclear matter at subsaturation density as a mat- 
ter has been subjected to a critical study for some time 
0, The aim is to understand the cv-clustering near 
the surface of heavy nuclei or the putative dilute alpha- 
condensate in light 4n nuclei. In astrophysical context, 
in following the evolution of the core-collapse supernovac, 
these studies have been extended to the case of warm nu- 
clear matter Q. The homogeneous low-density nuclear 
matter stabilizes as a mixture of nucleons and nucleon- 
clustcrs. It has lower free energy compared to that for 
nucleonic matter. The cluster composition is tempera- 
ture and density dependent, with increasing temperature 
or decreasing density, the population of heavier clusters 
tends to diminish leading to a mixture of nucleons and 
light clusters like d, t, 3 He and a 0,[1|. The properties of 
the clusterized matter undergo a major change, e.g., the 
incompressibility of clusterized nuclear matter is quite 
smaller compared to that for homogeneous nucleon mat- 
ter Q. This directly influences the collapse and bounce 
phase of the supernova matter. The symmetry energies of 
nuclear matter are also affected significantly when mat- 
ter gets clusterized @, H| ■ This has an important role in 
a better understanding of neutrino-driven energy trans- 
fer in supernova matter Q. The symmetry energy also 
influences the cluster composition in the crust of neutron 
stars and is thus instrumental in shaping the details of 
their mass, cooling and structure (loj . 

The equation of state (EOS) of warm dilute nuclear 
matter with only light clusters upto a has recently been 
investigated in the virial approach 0, El; inclusion of 
heavier clusters has also been made in the S'-matrix (SM) 
framework [l2|. These methods relate the calculations 
directly to the experimental observables like the bind- 
ing energies and the phase-shifts and thus, as such, are 
model-independent. They are usually taken as bench- 
mark calculations in the domain of low-density and high- 
temperature; they are understood to exhaust all the dy- 



namical information concerning the strong interactions in 
the medium. For an interacting quantum gas, the virial 
expansion, however, virtually ends at the second order. 
Formulation of higher order virial coefficients are very in- 
volved even at the formal level [l3[ , making it difficult to 
estimate the domain of validity of the virial series trun- 
cated at the second order. It may further be noted that 
the density should be dilute enough so that the concept 
of asymptotic wave functions as inherent in the virial ex- 
pansion should be meaningful. 

An alternate avenue could be to bypass the virial ex- 
pansion altogether and take recourse to nucleation in the 
framework of the mean-field model with a suitably chosen 
effective two-nucleon interaction that inherently takes an 
inclusive account of the scattering effects. Unlike the 
virial (S-matrix) approach which has direct contact with 
the experimental data, this method has indirect contact 
but it can be applied to relatively higher densities. With 
increasing density, a large number of different fragment 
species would, however, be formed that makes the nu- 
merical calculation very lengthy. Before attempting any 
full-blown calculation, it may then be worthwhile as a 
first step, to take only cv-clustering in the nuclear matter 
and to examine whether the model works in the low- 
density region where the bench-mark calculations exist. 
The present work aims towards that end. 

For the study of the so-mentioned cv-nucleon (aN) mat- 
ter, we have chosen the Thomas-Fermi prescription for 
the mean-field model and the finite range, momentum 
and density dependent modified Scylcr-Blanchard (SBM) 



effective interaction [141 ] . The properties we explore in- 
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elude the EOS of the aN matter, its symmetry energy, in- 
compressibility, a concentration etc. In Sec. II, the the- 
oretical framework for the mean-field and the S'-matrix 
approach is presented. Sec. Ill contains the results and 
discussions. Concluding remarks are given in Sec. IV. 



II. THEORETICAL FRAMEWORK 

Given an effective two-nucleon interaction, the prop- 
erties of the aN matter can be evaluated by exploiting 
the occupation functions of the n, p and alphas obtained 
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TABLE I: The parameters of the effective interaction (in 
MeV fm units) 



a 


Cu 


a 


b 


d 


K 


291.7 


910.6 


0.6199 


928.2 


0.879 


1/6 



system. In Sec. II A, some details of the effective interac- 
tion used are given. In Sec II B, theoretical formulation 
for obtaining the occupation functions from the Thomas- 
Fermi (TF) approximation is presented. In Sec. II C, 
expressions for various observables explored are given. 
In Sec. II D, a brief outline of the S'-matrix approach is 
made. 



A. The effective interaction 



The form of the SBM effective interaction v is 



v(r,p,p) = Ci. u [v 1 (r,p)+v 2 (r,p)}, 

r 2 



P 



V 2 



-(l-^)/(ri,r 2 ), 
d 2 [^i) + p(r 2 )] K /(ri,r 2 ), (1) 



with 



/(ri,r 2 ) 



e -|ri-r 2 |/a 

|ri - r 2 /a ' 



(2) 



The subscripts I and u to the interaction strength C re- 
fer to like-pair (nn or pp) and unlike-pair (np) interac- 
tions, respectively. The range of the effective interaction 
is given by a, b is the measure of the strength of the 
momentum dependence of the interaction. The relative 
separation of the interacting nucleons is r = ri — r 2 and 
the relative momentum is p = pi — p 2 ; d and k are the 
two parameters governing the strength of the density de- 
pendence and p(ri) and p(r 2 ) are the nucleon densities 
at the sites of the two interacting nucleons. The poten- 
tial parameters are given in Table I; the details for the 
determination of these parameters are given in [14 1. 
The incomprcssibility of symmetric nucleon matter 
is mostly governed by the parameter k; for the potential 
set we have chosen, the value of K ao =238 MeV. 

The equation of state of symmetric nuclear matter 
calculated with the SBM interaction is seen to agree 
extremely well fl5| with that obtained in a variational 
approach by Friedman and Pandharipande Q with 
t'14+TNI interaction. This interaction also reproduces 
quite well the binding energies, rms charge radii, charge 
distributions and the giant monopole resonance energies 
for a host of even-even nuclei ranging from 16 O to very 
heavy systems [l4[. Interactions of this type has been 
used with great success by Myers and Swiatecki [l|| in 



B. The occupation functions 

The self-consistent occupation probabilities of nucleons 
and alphas in qN matter at temperature T are obtained 
in the TF approximation by minimizing the thermody- 
namic potential of the system 



O = E - TS - VrN T - p a N a . 



(3) 



Here r represents the isospin index (n,p) . The quantities 
E,S, Pt i Pa i N T and N a are the total internal energy, en- 
tropy, nucleon chemical potentials, a chemical potential, 
free nucleon numbers and the number of a particles, re- 
spectively, in the system. Chemical equilibrium in the 
system ensures 



p a = 2(p n + p p ). 



(4) 



The total energy of the aN matter in TF approximation 
is 

E = J dr 1 dp 1 -^-n T (p 1 ) + i J dr 1 dp 1 dr 2 dp 2 

x[«!(|ri -r 2 |, |pi - p 2 |) + u 2 (|ri - r 2 |, 2p)][C l h T {p 2 ) 
+C u fi_ T (p 2 )]n T (pi) + i(C; + C u ) J dr 1 dp 1 dr 2 dp 2 

xn r ( Pl )n Q (p 2 )^ dr J dpfnf (pf )[«i(R', | Pi - (pf 

+p 2 )\)+v 2 (R\(^p T ,+p?))}\ 

f P 2 
+ / dridpi— — n Q (pi) 
J 2m Q 

+ (Ci +C U ) J dridpidr 2 dp 2 n Q (pi)n a (p 2 ) 



x 



dvdv' [ dp?dpfn?(pr)n?(Pi') 
x MR', |(pi + p?) - (P2 +pf'|) + v 2 (R', 2pf )] 



-N a B a 



(5) 



In Eq. (5), m T and m a are the nucleon and a masses, 
the first and fourth terms correspond to the kinetic en- 
ergy of the nucleons and alphas, the second and the fifth 
terms refer to the interaction energy among nucleons and 
among alphas, respectively and the third term is the in- 
teraction energy between nucleons and alphas. The var- 
ious space coordinates occurring in the third and fifth 
terms are shown in Figs. 1 and 2, respectively. 

These terms are evaluated in the single-folding and 
double-folding models. The last term is the binding 
energy contribution from the a particles. Here n T = 
-g^n T ,n a = pn a where n T and n a are the occupation 
probabilities for nucleons and alphas, respectively. Sim- 
ilarly, n"(p") = n"(pf) represents the occupation 
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FIG. 1: Space coordinates shown for nucleon (located at A) 
and alpha (with center at B) configuration. The origin of the 
coordinate system is at O and P is any arbitrary point within 
alpha. 




FIG. 2: Space coordinates shown for alpha-alpha configura- 
tion with O as the origin of the coordinate system. P and P' 
are arbitrary points within the alphas with A and A' as their 
centers. 



and pf is their intrinsic momentum inside the a. The 
space coordinates do not enter in the occupation func- 
tions n T and ?Tq CIS the system is infinite. For simplicity, 
the a particles are taken to be uniform nuclear drops 
with a sharp surface and hence the space coordinates do 
not also occur in hf. The notation J v refers to config- 
uration integral over the volume of a. The integral over 
p? is over the Fermi sphere of the nucleon momenta in- 
side the a particles. Since alphas are difficult to excite 
(the first excited state in a is ~ 20 MeV), they are taken 
to be in their ground states. All the other integrals are 
over the entire configuration or momentum space unless 
otherwise specified. It then follows that 



(n 



n a {p)dp = N a /V = p a , 
J hf(p)dp = 4/V a =pf, 



(6) 



where V is the volume of the aN system and V a = §7T-R^, 
with R a as the sharp-surface radius of the a drop taken 
to be 2.16 fm obtained from experimental rms charge- 
radius of a; pf is the density of the constituent nucleons 
of the a particles. The total baryon density pb is given 
by pi, = p + 4p a where p = ^Z T Pr is the density of the 
free nucleons and p a is the a-particle density. 
The total entropy of the aN system is 



s — 'S ' S T + S a , 

T 

where in the Landau quasi-particle approximation, 



(7) 



S T = 



J K(p) lnn r (p) 
+ (1 - n T (p)) ln(l - n T (p))]drdp, 



(8) 



and 



S a — 



1 



n a (p)ln n a (p) 



-(1 + n a (p))ln(l + n a (p))]drdp. 



(9) 



Minimization of fl with respect to n T and n a , remember- 
ing that 6n T (p) and 5n a (p) are separately arbitrary over 
the whole phase space, at the end yields 



dv- 



idp 2 {vi(| 



ri r 2 |, Pi p 2 



JL. 
2m T 

+v 2 (\r 1 - r 2 |,2p)|[Cin T (p 2 ) + C„n_ T (p 2 )] 
+Kd 2 (2p) K ' 1 Y J J dp , 1 dr 2 dp 2 

t' 

x.[C[h T ,(p 2 ) + C u h- T '(p2)]n T '(pi')f(T 1 ,r 2 ) 
+^(Q + C U ) J dr 2 dp 2 n a (p 2 ) 

x J drdp?n«( P n{«i(R',|Pi-(pr+P2)|) 
+v 2 (R' , (p + pf ) ) } + \ (Q + C u ) nd 2 (p + pf) ' 
J dp' 1 dp 2 h T ,(p' 1 )h a (p 2 )p? 

t' 

f f e ~ |R/|/a 
+T[lnn T ( Pl ) - ln(l - n T (pi))] - /i r = 0, 



(10) 



and 



n Jn)dn = NJV = n T 



Pi 



2(c, + a 



drodn->hr,(v>->) 
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x J drdprdT'dp? n?(p?)hf(p?) 

x{ Vl (R', |( Pl + pf) - ( P2 + p?')|) + v 2 (R',2pf )} 

+l(Q + C u )J2 J dpadpfftrfPaJfifCp?) 

x y^2^ dr{i> 1 (R',|p 2 -( Pl +p?)|) 
+U2(R',P + pr)}+T[lnn a (pi)-ln(l + n a (pi))] 



-(A» a + Sa) = 0. 



(11) 



Without any loss of generality, ri can be set equal to 
zero in Eqs. (10) and (11). The single-particle occupa- 
tion functions n T (p) and n a (p) for nucleons and alphas 
are determined from Eqs. (10) and (11), respectively. 
Eq. (10), after some algebraic manipulations can be writ- 
ten as 



A. 

2m T 



V»+p\V} + V; 



+T[lnn T (pi) - ln(l - n T ( P i))] - fx T = 0. (12) 

The momentum-dependent nucleon single-particle poten- 
tial V T (p) is given by 



V T (p) = V T °+p 2 V T 1 , 



(13) 



where V® is the momentum-independent part. Eq. (12) 
leads to 



n T (p) 



1 + exp 



/ 2 

/ p z 



V T ° + V 2 - p T IT 



\2m% 

where m* is the nucleon effective mass, 
— + 2K 1 



(15) 



and V 2 is the rearrangement potential coming from the 
density dependence of the interaction. Similarly Eq. (11) 
can be written as 



2m, 
-{pi a + B, 

which yields 
n a (p) = exp 
where 



K+pM + T[lnn Q (pi) - In(l + n a (pi))] 



^ + V0- { , a + Ba) \/T 



— + 

m a 



(16) 



(17) 



(18) 



is the a effective mass. V® is the momentum-independent 

j_ - r _ _ • l j.* l _ j j_* _ l t 7 / TrO i _2tt1\ 



the system. The nucleon and a masses are renormalized 
due to the momentum dependence in the interaction. 
The expressions for V® can be arrived at as, 

V T ° = -Ana 3 {1 - d 2 (2p) K } (C lPr + C u p_ T ) 
167r 2 a 3 " 



b 2 h 3 

X^3/2(?7-r) 

+d 2 (p + P TT 



d{2m* T Tfl 2 J 3/2 ( Vt ) + C u {2m*_ T Tf/ 2 



+ -I{Ci + C u ) PaP l 



<p 2 a > + < (pf) 2 > 



l> 2 



(19) 



The first two terms come from the interaction between 
free nucleons, the last term originates from the presence 
of alphas. In Eq. (19), I is the six-dimensional integral 
(see Fig. 1) 



f e -|r+R|/o 

/ = / dr dR- — — . 

Va J r + R/o 



(20) 



This integral can be evaluated analytically. The quantity 
< p 2 a > is the mean squared value of the a momentum 
in qN matter and < {pf) 2 > is the mean squared value 
of the constituent nucleon momentum inside the a. The 
value of < p 2 a > is 



<p 2 a > = (2m* a T)B 3/2 (r ]a )/B 1/2 (T la ) 
- 3m* T, 



and 



< (pf) 2 >c -(P£) 2 



(21) 



(22) 



where Pp is the value of the zero-temperature nucleon 
Fermi momentum inside a, taken to be 220.5 MeV/c, 
consistent with the a sharp surface radius. The Jkiv) 
and Bk(r]) are the Fermi and Bose integrals, 



Jk(v) 



x dx 
e (s-»7) + l ' 



and 



with 



B k (v) = I 
Jo 



x dx 



(23) 



(24) 



Vt = (p T - V r ° - V 2 )/T, 

Va = (p a + B a - V°)/T. (25) 

The expressions for , V 2 , V® and V* are given as 

Vt = ^[ClPr + C u p- T ] + hid + C u )^, (26) 



Vr 



Ana i Kd 2 {2pf~ 1 V[C„v + C u p- T ,]p T , 
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+ l m +c u )Kd 2 (p+p°r- i p?p aP , (27) 



V° = ^(Ci + C u )p? 2pfp a I a [d\2p?r - 1 



b 2 



+i[p{d 2 (p+ P ?r-i 



5 6 2 



and 



= \{Ci + C u )pf{2pf Pa I a + Ip}/b 2 . (29) 

In both and V 2 , the last term stems from the a- 
N interaction. The effective nucleon mass in pure nu- 
cleonic matter thus gets modified due to clusterization. 
The integral I a occurring in Eqs. (28) and (29) is a 
ninc-dimcnsional integral (see Fig. 2), 

f f f p-IR+r-r'l/a 

I a = I dr I dr' I dR,— — — — -— , (30) 



IR + r - r'|/a 



which can be evaluated numerically. If the alphas do 
not interpenetrate, the integral over R excludes the a 
volumes. 



and 



1 r 7T 

Pb 



+ ±(c l+ c v )i a pi(pf) 2 {d 2 (2p?r-i 

6m*T 6 (P$) 2 



b 2 



5 b 2 



(34) 



In the above equations, as stated earlier, pb (=p + 4p a ) 
corresponds to the total baryon density, p and p a are 
the free nucleon and a densities, respectively, in the aN 
system. 

ii) Entropy per baryon: The entropy per baryon Sb 
can be evaluated using Eqs. (8) and (9). It is additive 
and can be written as 



Sb = s N + s a , 



(35) 



where sn and s a are the contributions to entropy from 
free nucleons and alphas respectively. Their expressions 
reduce to 



Sat 



2 p5 



and 



Pa 
Pb 



(36) 



(37) 



C. Expressions for observables in TF 
approximation 

i) Energy per baryon : The energy per baryon e/, 
of the aN matter can be calculated from Eq. (5). It can 
be split into the following form, 



6fc — N + SaN + e c 



(31) 



Here cnn comes from the kinetic energy of the free nu- 
cleons and the interactions among them, e a N arises from 
the interaction among the free nucleons and the alphas 
and e aa stems from the kinetic energy of the alphas and 
the interaction among themselves. The expressions for 
them are 



= - VV [TJ 3 /2(Vr)/Jl/2(Vr){l - m*V?} 
Pb^ 



IK ], 



(32) 



— (Q + C u )I PaPa [[ 

i+d 2 (p+p?r}p 



(33) 



iii) Pressure of aN matter: Once the energy and 
entropy of the composite system are known, the pressure 
can be calculated from the Gibbs-Duhem thermodynamic 
identity, 



p = '^2 P r P T + P a P a ~ ft 



bPb, 



(38) 



where fb is the free energy per baryon, fb — eb — Tsb- 

iv) Incompressibility and the symmetry coeffi- 
cients: The incompressibility K can be computed from 
the derivative of pressure 



K = 9 



dP 
dp 



(39) 



The symmetry free energy and symmetry energy coeffi- 
cients Cf and Ce are calculated from 



Cf = - 



l(d 2 f b 



2\dX 2 Jx=o 



l(d 2 e b 
2 \dX 2 Jx=o 



(40) 



(41) 



where X is the neutron-proton asymmetry of the aN sys- 
tem. It is given as X = (p 7 b l — Pb)/Pb, where p\ l and p\ 

am i-]-tr\ fAf al nnnf rnn onrl rvrAf An rlnncifir rncnDptiirnlir 
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D. The S- matrix approach 

The relevant key elements of the S'-matrix framework 
17 1 as applied in the context of dilute nuclear matter 
8j, [12| are outlined in brief below. 

The grand partition function of an interacting infinite 
system of neutrons and protons can be written as 



Z= £ (tp) Z ((n) N Tr z 



n e 



-0H 



(42) 



Z.N=0 



where £ p 



and Cn 



arc the elementary fu- 



gacities with j3 = 1/T and /z's are the nucleonic chemical 
potentials. Here H is the total Hamiltonian of the sys- 
tem and the trace Ttz^n is taken over states of Z protons 
and TV neutrons. The partition function can be split into 
two types of terms [l7[ 



InZ = InZ 



(0) 
part 



In Z, 



(43) 



The first term on the right hand side corresponds to con- 
tributions from stable single-particle states of clusters 
of different sizes including free nucleons formed in the 
system; the second term refers to all possible scattering 
states. The superscript (0) indicates that the clusters be- 
have as an ideal quantum gas. In general, \nZ^J rt con- 
tains contributions from the ground states as well as the 
particle-stable excited states of all the clusters. The scat- 
tering term In Z scat may be written as a sum of scattering 
contributions from a set of channels, each set having to- 
tal proton number Z t and neutron number TV t . Since our 
interest in the present work is focused on aN matter, in 
lnZp^ rt , we include only the nucleons and the ground 
state of a; similarly in In Z scat , only the scattering chan- 
nels TV TV, aTV and aa are considered, so that 



In Z 



scat 



In Z NN + In Z aN + In 2 Q 



(44) 



Each of the terms in Eq. (44) can be expanded in the 
respective virial coefficients. Expansion upto the second- 
order coefficients are only considered. They are written 
as energy integrals of the relevant phase-shifts @, 0] . The 
partition function can then be written explicitly as 



In Z 



N 



+8C Q + 8b aa £ + 86 Q „Ca(C« + C P )]}, 



(45) 



where Xm= , h 



is the nucleon thermal wavelength, 

B a being the binding energy of a and 
/i Q = 2(// n + iip). The b nn , b np , etc.. are the temperature 
dependent virial coefficients @, Qjj]. The value of the 
virial coefficient b np has been adjusted so as to exclude 
the resonance formation of deuteron from n-p scattering 
to be consistent with our choice of the aN matter. 
The knowledge of the partition function allows all the 



given by 

P = T\nZ/V. (46) 
The number density pi is calculated from 

d lnZ\ 



Pi = d 



(47) 



where i stands for n,p, or a. Once the pressure, densities 
and chemical potentials are known, the free energy can be 
obtained from the Gibbs-Duhem relation. The entropy 
per baryon is calculated from 



Sb 



-(-) 



(48) 



which yields the energy per baryon as = fb + Tsb- The 
explicit expression for the entropy per baryon is 



T 



+Z£b' aa + 8U(n + ( P )b' an ] 



(49) 



The prime on the virial coefficients denotes their temper- 
ature derivatives. 



III. RESULTS AND DISCUSSIONS 

In the mean-field framework, the momentum 
and density-dependent finite-range modified Seyler- 
Blanchard force as scripted in Eqs. (1) and (2) has been 
chosen as the effective two-nucleon interaction in our 
calculations. To start with, we take baryon matter at 
a given density pb at a temperature T with an isospin 
asymmetry X. The unknowns are the free nucleon 
densities p n , p P and the a concentration in the matter. 
The three constraints are the conservation of the total 
baryon number, the total isospin and the condition of 
chemical equilibrium between the nucleons and alphas. 
Starting from a guess value for the a concentration, 
the unknowns are determined itcrativcly using the 
Newton-Raphson method. For our calculations, the 
masses of neutron and proton are taken to be the same, 
for a binding energy, the experimental value of 28.3 
MeV is used. For the evaluation of the aa-potential, 
the a-particles are assumed to be nuclear droplets with 
sharp boundary and that they do not interpenetrate. 

The calculations are done upto a baryon density 
/9b=0.01 fm~ 3 . To show the effect of temperature on 
different properties of the dilute matter, results are re- 
ported for temperatures T = 3, 5, and 10 MeV. In Fig. 3, 
the baryon fraction in a, Y a —Ap a / Pb (hereafter referred 
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0.002 0.004 0.006 0.008 
P h (fm 3 ) 



0.01 



FIG. 3: (color online) The a fraction Y a = 4p a /pb shown as 
a function of baryon density pb in TF and SM approaches at 
T=3, 5 and 10 MeV for symmetric matter (X=0.0) and asym- 
metric matter (X=0.2) in panels (a) and (b), respectively. 




0.002 0.004 0.006 0.008 0.01 
P b (fm 3 ) 

FIG. 5: (color online) Free energy per baryon F/A shown as 
a function of pb at T=3, 5 and 10 MeV in the TF framework 
for homogeneous nucleonic matter (blue lines) and aN matter 
(black lines). The red lines represent results from the SM 
approach. 




FIG. 4: (color online) The a fraction Y a displayed as a func- 
tion of asymmetry X at baryon density pb =0.001 (upper 
panel) and at 0.01 fm~ 3 (lower panel) at T= 3, 5 and 10 
MeV in TF and SM approaches. 



at the three temperatures mentioned arc shown for sym- 
metric (X=0) and asymmetric (X =0.2) nuclear matter 
in panels (a) and (b) , respectively. The black lines corre- 
spond to results obtained in the TF approximation [aN 
(TF)], the red lines refer to those in the SM approach 
[aN (SM)] with consideration of only n, p and a as the 
constituents of the baryonic matter. At low temperatures 
and higher densities, it is seen that alphas are the major 
constituents of the matter, with increasing temperature, 
the free nucleon fraction increases at the cost of a den- 
sity. At moderate asymmetry X=0.2, the a population 
is somewhat lower compared to that for symmetric nu- 
clear matter. In the temperature and density domain 
that we explore, the results from both the SM and TF 
approach are found to be quite close. The asymmetry 
dependence of a fraction Y a is displayed in Fig. 4 at two 
representative densities /9b=0.001 and 0.01 fm -3 at the 
three temperatures. With increasing asymmetry, the a 
concentration decreases, the decrease is more prominent 
at lower temperature. At the lower density (Fig. 4(a)), 
results for T=10 MeV are not shown as Y a is close to 
zero. 

In Fig. 5, the free energy per baryon for the homoge- 
neous nucleonic matter (denoted by N(TF)) and the aN 
matter in the TF approximation are presented in panels 
(a), (b), and (c) at T= 3, 5, and 10 MeV, respectively. 
The calculations presented refer to symmetric nuclear 
matter. The blue and black lines represent results for 
N(TF) and aN(TF). It is clearly seen that the cluster- 
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FIG. 6: (color online) Pressure P as a function of pt- The 
notations are the same as in Fig. 5. 



neous nucleonic matter. This is more prominent at lower 
temperatures, higher temperature tends to melt away the 
clusters. For comparison, results from the S'-matrix ap- 
proach are also presented. They are shown by the red 
lines, nearly indistinguishable from those from aN(TF). 
Fig. 6 displays the pressure of the baryonic matter. At 
lower temperatures (T=3 and 5 MeV), the nucleonic mat- 
ter shows the rise and fall of the pressure with density 
leading to unphysical region. For aN matter, however, no 
such unphysical region is observed in the density region 
we have studied. Both the TF and the SM approaches 
yield nearly the same value of pressure. At high temper- 
ature the a concentration becomes very less, the pressure 
in all the three approaches are then nearly the same in 
this density region. 

In Fig. 7, the effective masses of nucleon and a are 
shown as a function of density at the temperatures men- 
tioned. The nucleon effective mass is calculated for both 
nucleonic matter (blue line) and aN matter (full black 
line) in the TF approximation. The nucleon effective 
mass at a given pb in homogeneous nucleonic matter is 
always lower compared to that in clusterized matter. It 
is independent of temperature. In aN matter it nomi- 
nally decreases with temperature. At high temperature, 
the nucleon effective masses calculated in the homoge- 
neous and clusterized matter are nearly degenerate, with 
lowering of temperature, the degeneracy is lifted due to 
the increase in the a concentration. The effective a 
mass is shown by the dashed black lines. With increas- 
ing temperature, the medium effect on the a mass gets 
strikingly enhanced. This is due to the interplay of the 
temperature-dependent contributions from the aa inter- 
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FIG. 7: (color online) The nucleon (full black lines) and a 
(dashed black lines) effective masses shown as a function of pt 
at T=3, 5 and 10 MeV in the TF framework for aN matter. 
The blue lines refer to the corresponding nucleon effective 
masses for homogeneous nucleonic matter. 



and the second term within the braces in Eq. (29). 

The incompressibility of the baryonic matter as a func- 
tion of density is displayed in Fig. 8 at the three temper- 
atures. At very low density and higher temperature, the 
matter is mostly nucleonic in all the three approaches, 
so the incompressibility K is ~ 9T; this one sees at the 
lower densities considered at T =10 MeV in panel (c) of 
this figure. Even at this very high temperature, however, 
the nucleonic interactions have their role as the density 
increases; this results in the reduction of the incompress- 
ibility from the ideal gas value. At the lower tempera- 
tures (panels (a) and (b)), clusterization softens the mat- 
ter towards compression compared to homogeneous mat- 
ter (shown in the lower density region); increasing den- 
sity, however, pushes the homogeneous matter towards 
the unphysical region leading to negative incompressibil- 
ity. 

The symmetry energy coefficients Ce and Cf of the 
baryonic matter as a function of density are displayed in 
the left and right panels, respectively, of Fig. 9 at the 
three temperatures studied. The blue lines refer to cal- 
culations for the homogeneous matter, the black and red 
lines represent results for aN(TF) and aN(SM). Cluster- 
ized matter displays a marked increase in the symmetry 
coefficients noticed already earlier 0, |7j. The two ap- 
proaches to clusterization lead to the same values of the 
symmetry coefficients at lower densities, with increase in 
density the difference widens, more so at lower tempera- 
tures. 
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FIG. 8: (color online) The incompressibility K for baryonic 
matter shown as a function of pb at T=3, 5 and 10 MeV. The 
notations are the same as in Fig. 5. 



have been calculated with the assumption that the alphas 
do not overlap, they are mutually impenetrable spherical 
drops. This assumption relies on the fact that the alphas 
are very tightly bound and very hard to excite. To ex- 
plore the- effect of overlap in alphas, we consider a possi- 
bility of penetration with at best a 5 % overlap in volume 
(the value of I a in Eq. (30) then changes accordingly). 
Calculations have been repeated with this changed condi- 
tion. The so-calculated free energy per baryon, pressure 
and the a fraction Y a in the baryonic matter are pre- 
sented in panels (a), (b) and (c), respectively, of Fig. 10 
at T =3 MeV (the dot-dashed black lines) and compared 
with those calculated with the no-overlap condition (the 
full black lines) and also those from the aN(SM) ap- 
proach (the red lines). There is no significant change in 
the free energy or in a fraction, but the pressure changes 
perceptibly, particularly at higher density. The good 
agreement between the no-overlap aN(TF) calculations 
with those from the bench-mark aN(SM) shows the via- 
bility of the approximation of the impenetrability of the 
alphas. 



IV. CONCLUDING REMARKS 

Clusterization in warm dilute nuclear matter has been 
treated earlier in the virial approach or in the S'-matrix 
framework. These are model-independent parameter-free 
calculations. As explained in the introduction, these 
methods may have limitations at relatively high densi- 
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FIG. 9: (color online) The symmetry energy Ce (left panels) 
and symmetry free energy coefficients CV (right panels) shown 
as a function of at temperatures T=3, 5 and 10 MeV. The 
notations are the same as in Fig. 5. 
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FIG. 10: (color online) The free energy per particle, pressure 
and a fraction shown as a function of pb at T=3 MeV in panels 
(a), (b) and (c), respectively, for a drops with no overlap (full 
black lines) and with at best 5% overlap (dashed-dot black 
lines) in the TF approximation. The same observables are 
also shown in the SM approach (red lines). 
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ing with clusterized matter in a broadened density and 
temperature domain is suggested in the mean-held frame- 
work in the present paper. The suggested method may be 
lengthy at relatively higher densities where many differ- 
ent fragment species are formed, but it is straightforward. 
To explore its applicability in a wider domain, as a first 
step, we consider only n, p and a as the constituents 
of the matter at low densities and see how the results 
compare with those from the model-independent virial 
approach. 

We have chosen the SBM interaction that nicely re- 
produces the bulk properties of nuclear matter and of 
finite nuclei. We have calculated the a fraction, free en- 
ergy, pressure, incompressibility and the symmetry coef- 
ficients of this aN matter in this mean-field framework 
and find that all these results compare extremely well 



with those obtained from the iS-matrix method, partic- 
ularly in the low-density high-temperature regime. This 
gives one confidence in the applicability of this mean-field 
approach in dealing with the EOS of warm dilute bary- 
onic matter and the possibility of extending this method 
to higher densities. The price, however, is consideration 
of a larger number of fragment species and a numerically 
involved calculation. 
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